The Magnitude-Size Relation of Galaxies out to z ~ 1 



1,2 



Luc Simard^, David C Koo, S. M. Faber, Vicki L. Sarajedini, Nicole P. Vogt^'^, Andrew C 
Phillips, Karl Gebhardt, Garth D. Illingworth, K. L. Wu 

UCO/Lick Observatory, University of California, Santa Cruz, CA 95064, USA 

"a^": 
a^. 

^ : ABSTRACT 

xj ■ 

D ■ As part of the Deep Extragalactic Evolutionary Probe (DEEP) survey, a sample of 

1^ ■ 190 field galaxies (/814 < 23.5) in the "Groth Survey Strip" has been used to analyze 

G^ ! the magnitude- size relation over the range 0.1 < z < 1.1. The survey is statistically 

I I complete to this magnitude limit. All galaxies have photometric structural parameters, 

> ■ including bulge fractions (B/T), from Hubble Space Telescope images, and spectroscopic 

^ ■ redshifts from the Keck Telescope. The analysis includes a determination of the survey 

,__! ! selection function in the magnitude-size plane as a function of redshift, which mainly 

CN I drops faint galaxies at large distances. Our results suggest that selection effects play a 

qT ' very important role. 

G^ ■ A first analysis treats disk-dominated galaxies with B/T < 0.5. If selection effects are 

^ ! ignored, the mean disk surface brightness (averaged over all galaxies) increases by ~ 1.3 

Ph| mag from z = 0.1 to 0.9. However, most of this change is plausibly due to comparing low 

O ' luminosity galaxies in nearby redshift bins to high luminosity galaxies in distant bins. If 

+i . this effect is allowed for, no discernible evolution remains in the disk surface brightness 

j^ ! of bright {Mb < —19) disk-dominated galaxies. A second analysis treats all galaxies 

'^ I by substituting half-light radius for disk scale length, with similar conclusions. Indeed, 

at all redshifts, the bulk of galaxies is consistent with the magnitude- size envelope of 
local galaxies, i.e., with little or no evolution in surface brightness. In the two highest 
^ ', redshift bins {z > 0.7), a handful of luminous, high surface brightness galaxies appears 

that occupies a region of the magnitude-size plane rarely populated by local galaxies. 
Their wide range of colors and bulge fractions points to a variety of possible origins. 
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1. Introduction 

The field of galaxy evolution is rapidly changing 
under the wealth of new data coming from power- 
ful instruments such as the Hubble Space Telescope 
(HST) and 10-m class telescopes. Given the com- 
plexity of the problem, any attempt at understanding 
galaxy evolution must be deeply rooted in a quanti- 
tative and systematic study of fundamental galaxy 
structural parameters such as luminosity, size, mass 
and their associated scaling relations. High-redshift 
galaxies were first studied through their global lumi- 
nosity function. The landmark Canada-France Red- 



shift Survey (CFRS, [Lilly et al. 1995D showed that 
the luminosity function of blue galaxies evolves sig- 
nificantly over the redshift range 0.0 < z < 1.0 while 
that of red galaxies shows little or no evolution over 
the same redshift range. However, luminosity func- 
tions cannot paint a complete picture of the evolution 
of galaxies, as luminosity is only one of many galaxy 
structural parameters that may be evolving. 

Since individual galaxies cannot be observed at dif- 
ferent epochs, galaxy evolution can be operationally 
defined as the statistical change in galaxy properties 
with redshift. Correlations between observablcs such 
as total absolute magnitude M, half-light radius Rhi, 
disk scale length Rd, and rotation speed Vc are partic- 
ularly powerful diagnostics, as they involve key global 
properties of galaxies, namely luminosity, size, and 
mass. Furthermore, shifts in the zeropoints of such 
correlations are fairly insensitive to the local densi- 
ties of objects, and co-moving volume densities are 
therefore not required. 



So far, two gafaxy scaling relations have been ex- 



ness of ^ 1 mag at a given Rd by z — 0.6, both in the 
outer regions of clusters and in the field. Similarly, 
the effective surface brightness of early-type galaxies 
seemed to increase by 0.7 mag at z = 0.6. Using HST 
data, Roche et al. (1998) studied the surface bright- 
ness and size of 270 galaxies out to z = 4 taken from 
various redshift surveys. They found that the mean 
rest-frame blue-band surface brightness at a fixed ra- 
dius seemed to increase smoothly by 0.95 mag be- 
tween z = 0.2 and z = 0.9. They also concluded that 
a Pure Luminosity Evolution (PLE) model was not 
sufficient to explain the observed evolution in surface 
brightness and proposed a Size-Luminosity Evolution 
(SLE) model based on the evolution of different stellar 
populations as a function of galactoccntric distance. 
Lilly et al. (1998) concluded from ifST structural pa- 
rameters of 341 CFRS and LDSS galaxies that the 
increase in the disk surface brightness with redshift 
for galaxies with disk scale lengths greater than 4 
h^Q kpc was 0.8 ± 0.3 mag. This increase was not 
sufficient to account for the evolution of the overall 
luminosity function of galaxies seen in the CFRS sur- 
vey over the interval < z < 1. They suggested 
that smaller galaxies were responsible for most of the 
observed evolution, but they could not conclusively 
verify that hypothesis due to their magnitude and ra- 
dius limits. 

The question of whether the apparent evolution 
seen in galaxies at high redshifts is operating in both 
luminosity and/or size has received considerable at- 
tention in a recent flurry of theoretical activity aimed 
at explaining the formation and evolution of disk 
galaxies (Bouwens, Cayon & Silk 1997; Dalcanton 



Spergel fc Summers 1997|; [Mo, Mao fc White 1997 



plored at moderate redshifts. lire Al — 14 rcfation. 



also more commonly known as the 'iully-l'isher ('H') 



Mao, Mo & White 199?; Contardo, Steinmetz & 



Fritze-von Alvensleben 199^; Firmani fc Avila-Reese| 



relclion, has been studied to z ~ 1 by Vogt et al. 



(1996,1997). Based on only 16 galaxies, they found 
that: (1) massive, mature disks were present at z ^ 1, 
and (2) the change in the TF B-band zeropoint of 
these disks was at most 0.4 mag. Using ground- 
based CCD images, Schade et al. (1996a, b) presented 
Rd — Mdisk relations over the redshift range 0.1 < 
z < 0.6 for 351 late-type field and cluster galaxies 
and Re — Mb^ige relations for 166 early-type field and 
cluster galaxies over the same redshift range. Both 
samples were taken from the ground-based Canadian 
Network for Observational Cosmology (CNOC) sam- 
ple. The late-type galaxy magnitude-size relations 
showed an apparent increase in disk surface bright- 



1998|). Certain semi-ana lytical no-infall models (e.g.. 



Mo, Mao fc White 19971) predict that disk scale length 



should vary as (1 + z)^ in an fJ = 1 universe. Smaller 
disk scale lengths would be expected in the past due 
to the increase in the critical density Pcrit of the Uni- 
verse with look-back time. Dark matter halos are 
defined with a fixed overdensity {p'^ian ^ 200pcrit) 
in numerical simulations, and the radius of those ha- 
los must therefore decrease with increasing redshift 
to meet the overdensity criterion at high redshifts. 
Assuming that the mass and angular momentum of 
the gas settling into disks are fixed fractions of the 
mass and angular momentum of dark matter halos, 
smaller halo radii will translate into smaller disk scale 



lengths. Infall models (e.g., Chiappini, Matteucci & 
Gr ^tton 1997 ), on which the SLE model of Roche et 
al. (1998) is based, predict a milder disk size evolu- 
tion. Disk sizes at z ^ 1 in these models are typically 
only 20% smaller than at z = 0. Disk size is therefore 
important by itself to distinguish between these dif- 
ferent scenarios even though, as pointed out by Mao, 
Mo, & White (1998), it is likely that only kinematic 
information can decisively separate evolution in the 
mass and size of galaxies from evolution in their stel- 
lar mass-to-light ratio. 

One of the goals of the Deep Extragalactic Evolu- 
tionary Probe (DEEP; Koo 1998) survey is to com- 
pare the full 3D distribution of galaxy size, luminos- 
ity and rotat ion velocity at high r edshift to local re- 
lations (e.g., Burstcin ct al. 1997 ; hereafter BBFN) 
using a deep, homogenous sample of galaxies observed 
with the Hubble Space Telescope (HST) and the 10- 
m W.M. Keck Telescope. As a first step towards this 
goal, this paper focuses on the galaxy magnitude-size 
relation and seeks to extend previous work to fainter 
disk luminosities and higher redshifts. A major new 
focus is the role of selection effects. As this paper 



will show, selection ettects have had a major impact 



on the conclusions drawn from previous studies ot the 
magnitude-size relation. Further photometric data, 
internal velocities and full 3D structural distributions 
on larger samples of DEEP galaxies will follow in fu- 
ture papers. 

Section || describes the sample and observational 
measurements of galaxy sizes and magnitudes. Se- 
lection effects are quantified in Section 2.3. The 



magnitude-size relation as a function of redshift is 
presented in Section 0, and its implications for galaxy 
evolution are discussed in Section 0. Values of Hq = 
70 km s^^ Mpc^^ and go = 0.1 are used through- 
out this paper so that L* galaxies have / ^ 22.9 at 
z ~ 1 with 1" corresponding to 7.3 kpc. It should 
be noted that surface brightness is independent of 
go- AH magnitudes are on the Vega-based system 

{Isi4,Vega = hli.AB — 0.44). 

2. Observations and Analysis 

2.1. HST Structural Measurements 

The HST data are from two surveys, which we 
dub the "Groth Survey Strip" (GSS), taken under 
iJ5r programs GTO 5090 and GTO 5109. The GSS 
consists of 28 overlapping subfields taken with the 
HST Wide Field and Planetary Camera (WFPC2) 



and forms a "chevron strip" oriented NE to SW at 
roughly 1417-1-52 at Galactic latitude b ~ 60°. Each 
of 27 subfields has exposures of 2800 s in the broad 
V filter {F606W) and 4400 s in the broad / filter 
{F8UW). The 28th field is the Deep Field (J2000 
1417.5-^52.5), with total exposures of 24,400 s in 1/ 
and 25,200 s in /. The images were recalibrated "on- 
the-fly" through the Canadian Astronomy Data Cen- 
ter (CADC) standard pipeline and coadded with the 
IRAF^/STSDAS task CRREJ to remove cosmic rays. 

Object catalogs were produced with SExtractor 
version 1.0a (Berlin & Arnouts 1996) using a detec- 
tion threshold of l.ba and a minimum detection area 
of 10 contiguous pixels. SExtractor also produced a 
"segmentation" or "mask" image which was used to 
deblend galaxies from their nearby companions. Pix- 
els belonging to the same object are all assigned the 
same flag in this segmentation image, so the SExtrac- 
tor segmentation image was ideal for isolating objects 
in the surface brightness profile fits. 

The surface brightness profiles of galaxies in the 
DEEP/GSS sample were fitted w ith a PSF-con v olvcd 
2D bulge+disk mode l (GIM2D: [Simard 1998| |Mar- 
leau fc Simard 1998| ). The model had a maximum 
of twelve parameters: the flux Ftotai in data units 
(DU) integrated out to r = oo, the bulge fraction 
B/T (= Fbuige/ Ftotai), the bulge effective radius re, 
the bulge ellipticity e (e = 1 — b/a, b = semi-minor 
axis, a = semi-major axis), the bulge position angle 
of the major axis (fib, the exponential disk scale length 
Td, the disk inclination i (face-on = 0), the disk po- 
sition angle 4'd, the subpixel dx and dy offsets of the 
galaxy center, the background level b, and the Sersic 
index n. 

The first component (the "bulge" ) of the 2D sur- 
face brightness model was a Sersic profile of the form: 



E(r) = I]eexp{-fc[(r/re) 



l/n 



1]} 



(1) 



where S(r) is the surface brightness at r along the 
semi-major axis in linear flux units per unit area, and 
Eg is the effective surface brightness. The parame- 
ter k was set equal to 1.9992?i — 0.3271 so that r^ 
remained the projected radius enclosing half of the 
light in this component. The classical do Vaucouleurs 
profile has the special value n — A, and this value was 



^IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement 
with the National Science Foundation. 



chosen for the current analysis. This choice was mo- 
tivated by studies of bulge profiles in local galaxies. 
Locally, there is evidence that the bulges of late-type 
spiral galaxies may be better fitted by an n = 1 pro- 
file [^hereas^rig^|eliipticais^^ 



type spiral galaxies follow an n = 4 profile (|de Jong 
19^, pourteau et al. f 996| , [Andredakis f998| ). Local 
late- type galaxies with n — 1 bulges have B/T < O.f 



( |de| Jong f994 ). Since such bulges contain only fO% 
of the total galaxy light, low signal-to-noise measure- 
ments of late-type high-redshift galaxies make it very 
difficult to determine the Sersic index. On the other 
hand, n is more important for bulge-dominated galax- 
ies, and n = 4 is the expected value based on early- 
type local galaxies. The total flux in the Sersic bulge 
component is calculated by integrating Equation |l| 
from r = to infinity to obtain: 



Fbuige = 27rne*A:-'"r^r(2n)Ee 



(2) 



where F is the incomplete gamma function. For n = 
4, Fbuige = 7.2147rr2l]e. 

The second component (the "disk") was a simple 
exponential profile of the form: 



I](r) = Y.oexp{-r/rd), 



(3) 



where Eo is the central surface brightness. The total 
flux in the disk is given by: 



Fdisk = 27rr^Eo. 



(4) 



A PSF-deconvolved half-light radius r^i was also com- 
puted for each galaxy by integrating the sum of Equa- 
tions |I] and g out to infinity with the best fitting 
structural parameters. Even though the conventional 
"bulge/disk" nomenclature has been adopted in this 
paper, it should be kept in mind that this nomencla- 
ture does not say anything about the internal kine- 
matics of the components. The presence of a "disk" 
component does not necessarily imply the presence 
of an actual disk since many dynamically hot systems 
also have simple exponential profiles of the form given 
by Equation | ( |Lin fc Faber f983| , [Kormendy f985D . 
Likewise a "bulge" may represent a brightened center 
due to a starburst rather than a genuine dynamically 
hot spheroid. 

WFPC2 detector undersampling was taken into ac- 
count by generating the surface brightness model on 
an oversampled grid, convolving it with the appro- 
priate point spread function (PSF), shifting its cen- 



ter according to dx and dy, and rebinning the result 
to the detector resolution for direct comparison with 
the observed galaxy image. The PSF was generated 



by th e Space Telescope package TINY TIM ( |Krist 
f993| ) and subsampled to reproduce the pixel reso- 



lution of the Wide-Field camera. The best fitting 
parameter values were found with the Metropolis Al- 



gorithm (Metropolis et al. f953, 3aha & Williams 



f994), which Monte-Carlo samples parameter space 
to maximize the likelihood function. Best parame- 
ter values and their confidence intervals were deter- 
mined simultaneously by the algorithm as it refined 
its best solution. Figure |I| shows a 55" x 45" section 
of a typical WFPC2/F8f 4W GSS image before and 
after GfM2D processing. Figure || shows examples of 
GIM2D surface brightness fits for a set of five galaxies 
covering a wide range of bulge fractions and redshifts. 
Following the simulation procedure described in 
Section 3.4 of Marleau & Simard (f998), 400 simu- 
lated GSS galaxies were created with a wide range 
of structural parameters and analyzed in exactly the 
same way as real galaxies to test the reliability of 
the GSS parameter values measured with GfM2D. 
For simulated disk-dominated galaxies {B/T < 0.5) 
brighter than 1^14 = 23.5, the mean difference be- 
tween the measured and input disk magnitudes < 
AIsiA{disk) > was —0.02 with a{AIsiA{disk)) = 
0.20. For the same galaxies, the mean difference be- 
tween the measured and input log disk scale lengths 
< Alogrd > was 0.02, and a{Alogrd) was 0.05. Er- 
rors in bulge-l-disk fits are discussed in greater details 
in Marleau & Simard (1998). 

2.2. Keck/LRIS Spectroscopy 

GSS galaxies are currently being systematically 
surveyed by the DEEP team using the multi-object 
Low Resolution Imaging Spectrograph (LRfS, pke et 
al. 1995[) on the W.M. Keck If Telescope. DEEP/GSS 



galaxies were selected according to the magnitude cri- 
terion {V + I)/2 < 24. Spectra of a first group of 
231 galaxies were acquired between May 1995 and 
May 1997, and these are the galaxies used in this pa- 
per. Two gratings were used to cover a total spec- 
tral range of about 4500-9100 A depending on the 
exact position of the target on the mask: a 900 lines 
mm~^ grating (central wavelength ~ 5800 A, disper- 
sion = 0.85 A/pixel and resolution ~ 3-4 AfWHM), 
and a 600 lines mm"-'^ grating (central wavelength 
~ 7700 A, dispersion = 1.26 A/pixel and resolu- 
tion ~ 4-5 AFWHM). Typical total exposure time 



per target and per grating was 2700 seconds. Recti- 
fied, wavelength-calibrated, sky-subtracted 2D spec- 
tra were produced with a custom LRIS reduction 
package. One-dimensional spectra were optimally ex- 
tracted using the IRAF/APEXTRACT package. De- 
tails of the spectral data reduction will be given else- 
where. The final DEEP/GSS sample has a total of 
217 objects with both Keck redshifts and HST struc- 
tural parameters. The sample was cut down further 
to 190 galaxies by selecting only galaxies brighter 
than /814 = 23.5. The DEEP/GSS redshift sample 
is 100% statistically complete down to /8i4 — 23.5, 
meaning that redshifts were obtained for every object 
targeted spectroscopically above that limit. However, 
the redshift sample is not spatially complete since not 
all objects brighter than I^n — 23.5 were targeted. 
Furthermore, the basic photometric catalog may not 
be complete to /814 = 23.5 if low surface brightness 
galaxies are lost in the SExtractor selection process. 
This effect is discussed below but does not appear to 
affect our major conclusions. 

2.3. Selection Eflfects 

Since selection effects can conceivably mimic real 
evolutionary changes in the high-redshift galaxy pop- 
ulation, it is important to determine how they affect 
the DEEP/GSS sample. This section describes selec- 
tion effects both qualitatively and quantitatively. 

2.3.1. Observed Structural Parameter Distributions 

In order to determine whether the sample of galax- 
ies with redshifts is an unbiased representation of the 
photometric catalog population in the GSS, the ob- 
served structural parameter distributions of galaxies 
with Keck redshifts were first compared to those of 
725 general field galaxies measured in the same 6 fields 
of the GSS covered by the spectroscopic sample so 
far. These galaxies were photometrically selected us- 
ing SExtractor and then cut at /814 = 23.5, just like 
the sample of redshift galaxies. They should therefore 
be an equivalent sample, differing only in their lack of 
redshifts. Figure |^ compares the distributions of half- 
light radius, bulge fraction, magnitude, color and sur- 
face brightness for those galaxies in the photometric 
catalog population (6 GSS fields) to those with Keck 
redshifts. There are indeed no discernible differences 
between galaxies in the DEEP/GSS redshift sample 
and the photometric catalog population of the GSS 
fields. The bulge fraction distributions look nearly 
identical, and the color distributions also look similar. 



Finally, and most importantly, the surface brightness 
distributions show that the galaxies with DEEP/GSS 
Keck redshifts arc an unbiased sample of the surface 
brightnesses of the photometric catalog. A detailed 
analysis of the survey selection function is carried out 
in the next section. 

2.3.2. Determination of the DEEP/GSS Selection 
Function for Disk- Dominated Galaxies 

The observed distribution of galaxies in the Mbq — 
Rd plane as a function of redshift '^{Mbq, Rd, z) is 
the result of any inherent changes in the residentn 
galaxy distribution '^jj^Mbo , Rd, z) in that plane and 
of observational selection effects. Selection effects 
are likely to be significant given the wide range of 
sizes and surface brightnesses observed locally (p3en- 
dcr, Burstcin fc Faber 1992| , [Burstein et al. 1997|) . It 



is therefore important to carefully characterize selec- 
tion effects to disentangle them from real changes in 
■^u{MBo,Rd,z). The path from ^u{MBo,Rd,z) to 
*(Afs„, i?d, z) is given by: 



^{MBo,Rd, z) ^ Sps{MB,„Rd, z)Sup{MBo,Rd, z) 

^u{MBo,Rd,z), (5) 

where Mbo is the rest-frame B-band absolute magni- 
tude and Rd is the disk scale length in kpc (note the 
use of lower and upper cases to distinguish between 
apparent and intrinsic radii throughout this paper). 
The subscript UP stands for "Universe to Photomet- 
ric sample" , and the subscript PS stands for "Photo- 
metric sample to Spectroscopic sample" . The resident 
galaxy distribution 'i'ij{MB„, Rd, z) is not known a 
priori. Once the two selection functions in Equations 
have been characterized, their product (denoted Sus 
hereafter) shows the region of the Mbq — Rd plane 
where real galaxies would have been observed if they 
existed in that region at high redshift. 

The selection function Sup{MBo,Rd, z) contains 
the information needed to go from any sample of 
galaxies on the sky to the photometric catalog pro- 
duced with SExtractor and reflects the adopted SEx- 
tractor detection parameters (detection threshold in 



^It is very important to note the use of the term "resident" here 
and throughout the rest of the paper to refer to the intrin- 
sic galaxy population at a given redshift z. In the absence of 
real evolution in the galaxy population with redshift, all resi- 
dent populations would be the same as the local population of 
galaxies. 



sigmas, minimuni detection area, etc.). In practice, 
Sup{Mbo, Rd, z) is derived from the selection func- 
tion Sup{I?,iA,fd) determined as a function of the 
observed apparent magnitude /814 and the apparent 
disk scale length r^ in arcseconds. The transforma- 
tion Sup{hiA,rd) -^ Sup{MBQ-,Rd-,z) was made in 
each redshift bin using /c-corrections calculated with 
the median observed galaxy Vgoe — ^8i4 color of the 
B /T < 0.2 galaxies at that redshift and the assumed 
cosmology. 

Supilsi4, i~d) was constructed by generating 50,000 
B/T = galaxy models with structural parame- 
ter values uniformly covering the ranges: 16.0 < 
4i4 < 25.0, 0:'0 < ra < lO'.'O, 0° < i < 85°. 
Each model galaxy was added, one at a time, to 
an empty 20" x 20" section of a typical HST GSS 
image. "Empty" here means that no objects were 
detected by SExtractor in that sky section with the 
same detection parameters used to construct the ob- 
ject catalog. Using an empty section of the GSS en- 
sured that Sup{Isi4:, Td) was constructed with the real 
background noise that was seen by the detection algo- 
rithm. The background noise included read-out, sky 
and the brightness fluctuations of very faint galaxies 
below the detection threshold. This last contribution 
to the background noise is particularly hard to model 
theoretically, and the current approach bypassed this 
problem. SExtractor was run on each simulation with 
the same parameters that were used to build the SEx- 
tractor catalog. The function SupilsUi Td) was taken 
to be the fraction of galaxies successfully detected and 
measured by SExtractor at each value of (/8i4i ^d). 

The resulting S'[/p(/8i4, r^) is shown in the upper 
left-hand panel of Figure 0. Since the SExtractor de- 
tection algorithm depends mostly on surface bright- 
ness, it is easy to understand why Sup{l8i4:,Td) is 
bounded at low surface brightness by a line of nearly 
constant surface brightness, i.e., a straight line in the 
^Bo — Rd plane with slope close to —5. At the faint 
end, the boundary curves down as objects reach a 
total magnitude limit where they cannot meet the 
surface brightness detection threshold for any size. 
The upper center panel shows the magnitude distri- 
bution ^p(/8i4) of the general GSS population in the 
SExtractor catalog, and the upper right-hand panel 
shows the magnitude distribution ^s(/8i4) of galax- 
ies with Keck redshifts. The lower left-hand panel 
of Figure ^ shows the second observed selection func- 
tion, Spsilsu, Td) which is needed to go from the pho- 
tometric catalog ^_p(/8i4) to the spectroscopic sam- 



ple ^s(/8i4). Even though the DEEP/GSS sample 
is statistically complete down to /814 = 23.5 (i.e., a 
redshift was obtained for every galaxy targeted spec- 
troscopically) , Sps{I%i4,,Td) should not be expected 
to simply be a flat function down to that magnitude 
limit but will rather show a negative gradient since 
the fraction of galaxies in the photometric catalog 
targeted spectroscopically decreases with increasing 
magnitude. Sps{Is,iA, Td) is computed in the observed 
-^814 — I'd plane according to the rules described in 
Appendix A| and then transformed to the rest frame 
Sps{Mbo, Rd, z) in each redshift bin with the same k- 
corrections as for SupilsiAjTd)- Note that the cutoff 
at /814 — 23.5 has been applied. SpsihiA^rd) is flat 
at /814 < 19.0 since any galaxy above that magnitude 
limit would have been targeted spectroscopically. 

The final selection function Susilsu, fd), which is 
the product of Sup and Sps, is shown in the lower 
right-hand panel of Figure |j. The apparent sizes 
and magnitudes of galaxies with Keck redshifts above 
/814 = 23.5 are also plotted over the filled contours of 
Sus- Note that the observed density of points does 
not necessarily follow the contours of Sjjs since bright 
galaxies, though they have a greater probability of 
being detected, are also less numerous in the Uni- 
verse. The rest- frame counterparts of Susil&u, log 
rd) at each redshift are presented in the next section 
in conjunction with the magnitude-size relation for 
disk galaxies as a function of redshift. 

3. The Magnitude-Size Relation 
3.1. Disk galeLxies 

Figure H shows the relation between intrinsic disk 
scale length and disk rest-frame absolute B-band mag- 
nitude as a function of redshift for galaxies with 
B/T < 0.5 in the DEEP/GSS sample. The redshift 
ranges of the lowest three bins are nearly identical 
to those of Schade et al. (1996a,b). However, the 
Keck sample goes 1.5-2.0 magnitudes deeper than 
their CNOC sample and also provides two bins at 
higher redshifts between z = 0.70 and z = 1.10. The 
error bars on the points are the 99% confidence in- 
tervals of the /814 and rd measurements. The open 
triangles are objects with structural parameters mea- 
sured with the same software as the rest of the objects 
but whose redshifts come from the CFRS survey. The 
redshift range covered by all five bins corresponds to 
a total lookback time interval of 5.8 Gyrs {qq = 0.1). 

The long-dashed line is the Freeman relation for an 



exponential disk B-band central surface brightness of 
21.65 mag arcsec"^ (Freeman 1970). The disk sizes 
and luminosities of most galaxies in the two lowest 
redshift bins cluster near this relation. Galaxies in 
the third redshift bin also cluster around the Freeman 
relation, although a substantial number of objects be- 
gin to exhibit systematically higher central surface 
brightness than the Freeman value. Major deviations 
from the Freeman relation are seen above z — 0.70, 
with evidence for a class of very small, high-surface 
brightness objects at z > 0.90. It is interesting to 
observe that a significant number of galaxies remain 
on the Freeman locus even at z > 0.90, consistent 
with the existence of a population of Freeman sur- 
face brightness disks even at this early epoch. The 
large spread in surface brightness at z > 0.90 may be 
due either to a strong differential disk evolution that 
leaves some galaxies unaffected or to a new galaxy 
population making its entrance at z > 0.90. Both 
possibilities are discussed further in Section ^. 

The magnitude-size relations shown in Figure |5| 
cannot be properly interpreted without taking selec- 
tion effects into account. In Figure g, the rest-frame 
selection function Sjjs for each redshift bin is shown 
as the shaded contours, and the same galaxies as in 
Figure ^ are replotted to compare their location to 
the predictions of the selection function. The Free- 
man disk relation is also replotted for reference. Two 
important observations can be made from Figure o. 
First, very few disks are observed at low redshifts 
in the ranges —20 < MBo{disk) < —18 and 0.5 < 
log Rd < 1.2 (significantly above the Freeman rela- 
tion) even though they are detectable in principle. 
Second, the envelope of the selection function in the 
two highest redshift bins comes close to the Freeman 
locus, and this situation has a very serious implica- 
tion: galaxies whose surface brightness may have been 
fainter at those redshifts are not detectable in the 
DEEP/GSS sample. Hence, this sample cannot be 
used to distinguish between a scenario in which mean 
surface brightness is higher at high redshift and one 
in which mean surface brightness remains unchanged 
but the spread in surface brightness about that mean 
value increases. 

Although the present sample may be somewhat bi- 
ased to high surface brightness at high redshift, this 
effect would appear to be even more pronounced in 
the Schade et al. (1996a,b) and Roche et al. (1998) 
samples. This is shown in Figure |^, which illustrates 
loci in the magnitude-size plane occupied by the vari- 



ous samples as a function of redshift. The DEEP/GSS 
loci are taken from Figure H and are shown by the 
solid boxes. Loci for the other samples have been 
drawn by eye from analogous plots in the published 
papers. These boxes cannot substitute for a quanti- 
tative treatment, but even this rough comparison is 
instructive. 

Systematically, it is seen that the loci of the Schade 
et al. and Roche et al. samples lie at higher surface 
brightness relative to the DEEP/GSS sample, and 
that this offset increases with redshift. Put differ- 
ently, the first two samples fail to include many low 
surface brightness objects of the sort actually detected 
by DEEP/GSS. Quantitatively, at z ~ 0.7, the up- 
per boundary difference between the first two samples 
and DEEP/GSS is about 0.3 dex in log Rd- Loss of 
these galaxies would translate to a shift of about 0.75 
mag in mean surface brightness if the DEEP/GSS box 
were uniformly populated vertically - about the same 
magnitude as the surface brightness evolution actu- 
ally found in the other two samples. Figure thus 
raises the possibility that the surface brightness evo- 
lution seen by Schade et al. and Roche et al. might 
have been caused in large part by the loss of low sur- 
face brightness objects at high z. 

Figure M also shows boxes marking the distant com- 
pact galaxy sample of Phillips et al. (1997). This 
sample represents the top 35% of galaxies in surface 
brightness to Isu — 23.5 around the Hubble Deep 
Field. The location of the Phillips et al. sample in 
Figure M is consistent with its compact high surface 
brightness nature. Further comparison is made to this 
sample below. 

As a first step to study in greater detail whether 
the magnitude-size relations shown in Figure ra ex- 
hibit real evolution or not, surface brightness distri- 
butions were calculated for each redshift bin. The 
disk rest-frame B-band central surface brightness of 
each galaxy was computed using the relation: 



Hob = Mbo + 5.01ogi?d + 38.57, 



(6) 



where Mbq is the disk rest-frame B-band absolute 
magnitude, and Rd is the disk scale length in kpc. 
Equations is valid for pure exponential disks, and the 
units of fiQg are mag arcsec"^. The resultant surface 
brightness distributions shown in Figure are raw 
without correction for selection effects, and the area 
under each distribution is normalized to unity. The 
mean of each distribution is shown as a large filled 
circle, and the error bar attached to each circle is 3 



times the standard error of the mean. The vertical 
dotted hnes mark the canonical Freeman disk central 
surface brightness of 21.65 mag arcsec"^. There is a 
systematic increase in the mean disk central surface 
brightness from /zo^ = 21.7 mag arcsec"^ at z = 0.20 
to /ios = 20.4 mag arcsec"^ at 2; = 0.80. This raw 
increase of 1.3 mag is similar to the effect seen by 
Schade et al. and Roche et al. over the same redshift 
range. 

Next, the surface brightness distributions were re- 
calculated by taking the survey selection function into 
account. Let c/jij {fJ-Og. . )dfioB. - be the number of galax- 
ies in the j*^ bin of the surface brightness distribution 
for the i*'* redshift bin. Adjusted surface brightness 
distributions in Figure were calculated using the 
equation: 



</>Jj(/^0B,,)'^M0i5, 



Sus{Rd,,„M,jk, 0-90 <z) 



Z^ 



Sus{Rd,,k,Mijk,z.i) 



(7) 
where the index k runs over all the galaxies in the i*'* 
redshift bin with surface brightness between /xq^. . and 
/Uob + rf/J-Ofi. - • Each object was thus de- weighted by 
the inverse ratio of the value of the selection function 
in its redshift bin and magnitude-size position to the 
value of the selection function in the highest redshift 
bin at the same magnitude-size position. This weight- 
ing scheme effectively applies the selection function of 
the highest redshift bin uniformly to all the redshift 
bins. The area under each distribution in Figure H 
was again normalized to unity. The parameters of 
the surface brightness distributions with and without 
selection function corrections are given in Table ^. 
These distributions show that there is no detectable 
systematic change in the mean disk central surface 
brightness over the redshift range 0.1 — 1.1 when a uni- 
form selection function is applied to all bins. Why has 
this occurred? Previous discussion hinted that loss of 
low surface brightness galaxies at high z might be im- 
portant. However, the major consequence of apply- 
ing the highest redshift selection function to all bins 
is actually the deletion of intrinsically faint galax- 
ies below Mbq = — 19 mag from the nearer redshift 
bins. Inspection of Figure shows that the remain- 
ing, brighter galaxies tend to lie below the Freeman 
relation, i.e., at higher surface brightness, at all red- 
shifts. In other words, once comparison is restricted 
to galaxies of similar (i.e., bright) absolute magni- 
tude, no significant evidence for evolution in mean 



disk surface brightness remains, at least for galaxies 
in the range MBo{disk) = — 19 to —22 mag. 

3.2. Global Population 

The analysis so far has been restricted to the disks 
of disk-dominated galaxies with B/T < 0.5. To ex- 
tend the analysis to the whole distant galaxy popu- 
lation, a local magnitude-size relation must be cho- 
sen that covers as many local galaxies as possible. 
The choice of a local magnitude-size relation is clearly 
key to interpreting the relation at higher redshifts. 
Schade et al. ( 1996a, b) used fits to galaxies in the 
cluster Abell 2256 (z — 0.06) as a local reference. 
Roche et al. (1998) used a set of four galaxy type- 
dependent relations (Freeman, E/SO, Sab-Sbc, and 
Scd-Im). While these relations are useful, they do 
not represent very well the wide range of local galaxy 
photometric parameters over all Hubble types. 

To broaden the approach, half-light radii Rhi are 
introduced as a measure of galaxy size. Half-light 
radii are available for all local Hubble types in quan- 
tity in the RC3 catalog. Specifically, the sample of 957 
local galaxies extracted from the RC3 by Burstein et 
al. (1997; BBFN) is used here. Figure |i^ shows half- 
light radius versus Mbq for this sample. The data are 
divided into four broad morphological classes: ellip- 
ticals (E), early- type spirals (Sa-Sbc), late- type spi- 
rals (Sc-Sdm), and irregular galaxies (Sm-Irr). The 
Freeman relation is again shown as a reference (now 
expressed using half-light radius). This large dataset 
demonstrates how significant the total spread in the 
local magnitude-size relations is: ~ 1.5 mag for the 
bright ellipticals and ~ 3.0 mag for the late-type spi- 
rals. More significant is the fact that the different lo- 
cal Hubble types differ systematically in radius (and, 
therefore, in surface brightness) at a given absolute 
magnitude. Figure |l^ shows that, at Mbo — —22 
mag, bright Sa-Sc spirals are nearly 2 magnitudes 
brighter in mean surface brightness than ellipticals; 
however, at Mbq = — 18 mag, they are 2 magnitudes 
dimmer; Sm and Irr galaxies at Mbq = —18 are a 
further 0.5 magnitude dimmer yet than Sc's. 

These systematic effects are of crucial importance 
in predicting what distant surveys will show. Sup- 
pose, for the sake of argument, that the distant galaxy 
population is the same as that seen locally. The dis- 
tant selection function will then cut through the man- 
ifold of galaxies in different locations at different red- 
shifts, turning up galaxies in proportion to their res- 
ident number densities at each redshift. At nearby 



redshifts, the selection function will turn up dim but 
abundant Sm and Irr types. At distant redshifts, it 
will turn up a mixture of rarer but brighter E's and 
early-type spirals. Mean surface brightness averaged 
over all detected galaxies will thus tend to "walk" 
from low to high values with distance, precisely as 
seen in the raw data in Figure @. (It may even walk 
at fixed absolute magnitude, depending on the rel- 
ative volume densities of different Hubble types and 
how their surface brightnesses differ at that particular 
magnitude.) 

Finally, it should be stressed that the BBFN cat- 
alog in Figure M is not a volume-limited sample, nor 
does it fairly sample the full range in surface bright- 
nesses seen locally, from high surface brightness star- 
bursting HII and compact narrow emission-line galax- 
ies on the one hand, to low surface brightness dwarf 
spheroidals on the other. Comparisons must therefore 
be made with caution, taking note of these deficien- 
cies. These issues are discussed further below. 

Measured half-hght radii of all DEEP/GSS galax- 
ies are plotted against total rest-frame B-band abso- 
lute magnitudes in Figure |l^. Filled circles are disk- 
dominated galaxies {B/T < 0.5), and open circles are 
bulge-dominated galaxies. The solid line is the Free- 
man relation. Figure O again shows that a significant 
number of galaxies remain on or close to the Freeman 
relation at all redshifts, and further that such galaxies 
are largely disk-dominated. Bulge-dominated galax- 
ies appear to lie at higher surface brightnesses than 
disk-dominated galaxies, especially at the highest red- 
shifts. 

To facilitate further comparison. Figure n^ directly 
overplots the distant half-light radii from Figure O 
on top of the total local catalog from Figure |l^ (tiny 
dots). As expected, at low redshift, the DEEP/GSS 
sample is dominated by intrinsically faint but abun- 
dant late-type galaxies deep in the luminosity func- 
tion. These galaxies are largely disk-dominated, with 
B/T < 0.5 (solid dots. Figure |Il|) and occupy the 
same magnitude-size locus as local Sm-Irr galaxies 
(Figure [1C|) , which also hav e pure exponential profiles 
( Gallagher fc Hunter 1984| ). The excellent agreement 
between the two loci is not surprising given the small 
difference in epochs between the two samples. 

As redshift increases, the observed magnitude-size 
envelope begins to shift, both in surface brightness 
and also in absolute magnitude. However, out to z = 
0.70, the projected distributions in the magnitude- 
size plane remain roughly consistent with those pre- 



dicted by the local envelope. There may be a slight 
shift to higher surface brightness at Mb,, ^ —20 in the 
z = 0.50-0.70 bin, but the number of points is small. 
Overall, the impression is one of no great change from 
local galaxies out to z = 0.70. 

Beyond z = 0.70, there is a hint of change. First, 
the range of surface brightnesses in the z — 0.70-0.90 
bin appears to be larger than that seen locally, and a 
handful of very high surface brightness objects (HSBs 
hereafter) has appeared. The latter effect cannot be 
due to observational errors on the radii, as they are 
quite small (cf . Figure 0) . Denoted by the open cir- 
cles in Figure ll^, these HSBs increase still further in 
number in the z — 0.90-1.10 bin, where objects dis- 
tinctly more compact than the local BBFN sample 
are seen. Remarkably, HSBs span the whole range 
of bulge fractions and rest- frame colors and therefore 
may not be restricted to a single structural type. 

HSBs aside, the bulk of galaxies at z > 0.7 are no- 
table more for their resemblance to the local sample 
than for their differences. This statement is tentative, 
as no quantitative comparison has yet been made be- 
tween the near and distant samples taking resident 
local densities of the various structural types properly 
into account. It is also conceivable that further differ- 
ences will emerge once velocity widths are measured 
(a survey of these is in progress). For now, the present 
results suggest simply that no great global evolution 
has occurred in the surface brightnesses of luminous 
galaxies since z ^ 1. This conclusion is based on a 
considerably deeper and more complete sample than 
heretofore available, as Figure M shows. 

Strictly speaking, one should also carry out selec- 
tion function experiments on half-light radii of the 
various Hubble types, analogous to those carried out 
for disk radii in Figure ra. However, the main ef- 
fect of these experiments for bright galaxies was to 
limit the sample in magnitude, not in surface bright- 
ness. Thus, neither the relative constancy in surface 
brightness nor the appearance of luminous high sur- 
face brightness galaxies are likely to stem from as-yet- 
uncalibrated selection effects. 

4. Discussion 

The inclusion of selection effects in the analysis of 
the disk galaxy magnitude-size relation at high red- 
shift has led to a number of important results. 

First, even though the DEEP/GSS survey goes 
nearly two magnitudes deeper than previous studies. 



the faint envelope of the survey selection function in 
surface brightness approaches the local Freeman rela- 
tion at redshifts greater than z — 0.70. This means 
that the magnitude-size relation of past and present 
samples cannot be used by itself to favor a scenario in 
which the mean surface brightness of disks brightens 
systematically with redshift over a scenario in which 
the spread in disk surface brightness increases. Al- 
though theoretical expectations for the characteris- 
tic surface brightness of "normal" bright disk galax- 
ies are comparable to the canonical central Freeman 
value of 21.7 ± 0.3 mag arcsec"^, the spread about 
that value could be as large as ± 3.4 mag depending 
on the range of disk masses and mass -to-light ratios 
( Dalcanton, Spcrgcl fc Summers 1997 ). The disks on 
the faint end of that surface brightness distribution 
would quickly be excluded from the DEEP/GSS and 
other samples at relatively low redshifts. 

Second, selection effects produce an apparent sys- 
tematic increase in disk mean surface brightness that, 
at first sight, looks like luminosity and/or size evolu- 
tio n (Figure p|). However, there is no such systematic 



inc 'ease below z — 1.1 when the selection function 



is used to weigh the disk surface brightness distribu- 
tions as a function of redshift. According to the se- 
lection function Syg , disks with absolute magnitudes 
brighter than Mb^ = — 19 remain detectable over all 
redshift bins, and these are the disks that effectively 
contribute to the weighted surface brightness distribu- 
tions. Significant evolution in the surface brightness 
of such disks is not apparent, though the numbers 
are still small. Selection effects are too important for 
fainter disks to determine whether or not their surface 
brightness is evolving with redshift. 

Third, it is apparent from Figures o and Rfi that 
a significant number of galaxies remain close to the 
canonical Freeman relation at all redshifts. This pop- 
ulation may be the same as that detected in the TuUy- 
Fisher zeropoint studies of distant galaxies by Vogt et 
al. (1996,1997). The Vogt galaxies were selected to 
include the largest (i.e., best resolved) morphologi- 
cally normal spirals at each apparent magnitude, a 
criterion that favors lower surface brightness galax- 
ies in each redshift bin. It is probable that the Vogt 
et al. studies also targeted late-type spirals prefer- 
entially, since they have stronger emission lines. The 
constancy of the TF zeropoint for these particular 
objects would then not be surprising, since the star- 
formation timescales of late-type galaxies are quite 
long {tsfr > 7 Gyrs, Bruzual fc Chariot 1993 ), and 



mass-to-light ratios for such stellar populations evolve 
only slightly. 

These long timescales are a feature of the Pure 
Luminosity Evolution model of Roche et al. (1998), 
which predicts that the surface brightness of late- type 
spirals should remain constant out to z = 1. On the 
other hand, their Size-Luminosity Evolution model 
predicts that the same late-type galaxies should in- 
crease nearly 1.5 magnitudes in surface brightness to 
the same redshift (Figure 7, Roche et al. 199S). This 
latter model does not appear to be consistent with 
the DEEP/GSS data. Indeed, according to Figure |ll], 
such a large brightening, if present in all disk galaxies, 
would considerably decrease the number of galaxies 
on the Freeman relation. This deficiency of galaxies 
on the Freeman relation could not be compensated 
by increasing the surface brightness of any resident 
Irr population, even if present. The internal kinemat- 
ics of Irr galaxies cannot reproduce the observations 
of Vogt et al. The most massive Irr galaxies have ro- 
tation velocities lower than 100 km s~^, and typical 
Irr galaxies h ave Vc ~ 50 — 70 km s^^ ( [Gallagher & 



Hunter 1984) 



As noted earlier, the last redshift bin (z > 0.9) 
contains about nine luminous, high surface bright- 
ness (HSB) objects (shown as open circles in Fig- 
ure |lj) that occupy a region of the magnitude-size 
plane rarely populated by local galaxies. The ex- 
act number of these aberrant objects depends on 
the boundary adopted for the envelope of normal lo- 
cal galaxies; a conservative, high surface brightness 
boundary was adopted here, minimizing the number 
of HSBs. The nature of the HSBs will be the subject 
of a future paper. For now, it is sufficient to note that 
the present HSBs occupy essentially the same region 
of the magnitude-size plane as the distant compact 
sample of Phillips et al. (1997), as shown in Figure R. 
Like the present HSBs, that sample exhibited a mix- 
ture of colors and emission-line strengths, and many 
also had structural parameters similar to the distant 
compact narrow emission-line galaxies of Guzman et 
al. (1998). Based on colors and spectra, Phillips et al. 
concluded that their population of HSBs was rather 
heterogeneous, though with a tendency to low internal 
velocity widths compared to local galaxies of similar 
absolute magnitude. The present data give further 
support to the notion that HSBs are a mixed popu- 
lation. Furthermore, it is also worth noting that the 
magnitude-size locus of the z ~ 3 galaxies observed 
by Lowenthal et al. (1997) overlaps with some HSBs 
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(see dashed box in Figure 12). Although this overlap 
raises interesting possibilities, it is by itself insuffi- 
cient to directly link some HSBs to z ~ 3 galaxies. 
Potential links between these two types of galaxies 
will be explored in a future paper using a full set of 
structural parameters and emission linewidths. 

At least three explanations for HSBs suggest them- 
selves: (1) They may represent an entirely new pop- 
ulation of galaxies not found locally and making its 
appearance at z '^ 1. (2) The local magnitude-size 
dataset used as a comparison here may be incom- 
plete. Indeed, BBFN warn that rare types of local 
HSB galaxies such as starburst HII galaxies may be 
missing from their sample. The blue colors of some 
HSBs strongly suggest that at least some are star- 
bursts. (3) HSB galaxies may have had high surface 
brightness in the past but have faded to become a nor- 
mal population locally. For example, at a fixed radius, 
most HSB galaxies would have to fade by at least 1 
mag to fall within the local magnitude-size envelope 
today. This is consistent with the passive evolution of 
a very early-type population, and it would nicely ex- 
plain red HSBs. In view of the heterogeneous nature 
of HSBs, all three scenarios may be at work. 

5. Conclusions 

The HST structural parameters of 190 galaxies 
ihii < 23.5) in the Groth Survey Strip with DEEP 
Keck redshifts have been used to build a magnitude- 
size relation in five redshift bins over the range 0.10 
< z < 1.10. Selection functions were quantitatively 
determined for all five redshift bins. The major con- 
clusions are: 

1. If selection effects are ignored, the mean surface 
brightness of disks in the sample appears to increase 
systematically by ^ 1.3 mag from z = 0.1 to z = 
0.9. Such an increase was also found by Schade et al. 
(1996b) and Roche et al. (1998). 

2. However, if the survey selection functions are 
used to correct the surface brightness distributions in 
the different redshift bins, then there is no detectable 
change out to z = 0.9 in the mean surface bright- 



ness of disks brighter than Mbq 



-19. These disks 



are the faintest objects not significantly affected by 
selection effects in the highest redshift bin. 

3. A number of disks remain on the canonical 
Freeman relation in all redshift bins. This is con- 
sistent with the lack of evolution seen in the zero- 



point of the high-redshift Tully-Fisher relation of Vogt 
et al., and it does not support a 1.5 mag brighten- 
ing in surface brightness predicted for all late-type 
galaxies by the Size-Luminosity Evolution model of 
Roche et al. (1998). On the other hand, the constant 
surface brightness of late-type galaxies predicted by 
their Pure Luminosity Evolution model agrees with 
the DEEP/GSS observations of these disks. 

4. At all redshifts, two important factors conspire 
to produce an apparent increase in global galaxy sur- 
face brightness out to z = 0.9: the relative volume 
densities of different Hubble types and the survey se- 
lection function. Lower surface brightness late-type 
and irregular galaxies will dominate lower redshift 
bins due to their higher volume densities. Higher sur- 
face brightness early-type galaxies will dominate the 
highest redshift bins as a result of increasing survey 
volume and the selection function cutting different 
slices through the resident galaxy surface brightness 
distribution at that redshift. So, for example, even 
if the distant galaxies were identical to those seen 
locally, average galaxy surface brightness would still 
shift to higher values at higher redshifts, in qualita- 
tive agreement with what is seen. Systematic shifts 
in galaxy surface brightness versus redshift thus may 
not necessarily represent real evolution in the radii or 
luminosities of the resident galaxy population. 

5. Nine galaxies at redshifts z > 0.9 occupy a re- 
gion of the magnitude-size plane rarely populated by 
local galaxies. As will be shown in a forthcoming pa- 
per, these luminous, high surface brightness galaxies 
exhibit a wide range of rest-frame colors and bulge 
fractions. This suggests that a variety of scenarios 
might be required to explain their origin: a new pop- 
ulation at z ~ 1 not found locally, a rare population 
missed by local surveys, and/or a familiar local pop- 
ulation that has faded by at least 1 mag from z '--^ 1 
to the present epoch. 

Further work is needed to fully understand the 
magnitude-size relation of the global galaxy popula- 
tion out to z '^ 1. Prime requirements are a truly 
volume-limited sample of local galaxies plus comov- 
ing volume densities for the distant samples. The 
size of the distant sample must also be increased to 
study the statistical properties of the HSB population 
and determine its origin. Finally, linewidths must be 
obtained to provide full structural information, in- 
cluding masses. These additional parts of the global 
structural puzzle will be tackled in future papers. 
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Fig. 1. — Upper Panel: 55" x 45" section of a typical 
WFPC2/F814W GSS image before GIM2D process- 
ing. Exposure time was 4400s. Galaxies in this image 
section span a wide range of bulge fractions. Lower 
Panel: GIM2D residual galaxy image of the same im- 
age section as shown above. GIM2D produced this 
residual image by subtracting the best bulge-|-disk 
surface brightness model from each galaxy detected 
by SExtractor. 

Fig. 2. — Examples of GIM2D surface brightness fits 
for a set of five galaxies covering a range of bulge frac- 
tions and redshifts. The left-hand panels are the orig- 
inal F814W galaxy postage stamp images analyzed 
by GIM2D, the central panels are the PSF-convolved 
best fitting surface brightness models, and the right- 
hand panels are the residual images. Redshift and 
measured bulge fraction are given for each object. 
These five galaxies have 19.9 < Isu < 22.3. 

Fig. 3. — Observed structural parameter distributions 
for galaxies with DEEP Keck redshifts {solid line) and 
for an identically selected sample of galaxies in the 
general GSS population {dotted line). From left to 
right and from top to bottom: PSF-deconvolved half- 
light radius rhi, observed FSIAW bulge fraction B/T, 
FSMW total magnitude Isi4, total Vgoe —-^814 galaxy 
color, observed F814M^ effective surface brightness 
/Xesi4 in mag arcsec"^. Magnitudes and colors are in 
the Vega zeropoint system. The spectroscopic sam- 
ple shares the same characteristics as the general field 
population. 

Fig. 4. — From left to right and from top to bottom: 
SExtractor selection function Sup{Isi4:,i"d)i magni- 
tude distribution ^p(/8i4) of galaxies in the photo- 
metric catalog, magnitude distribution ^s(/8i4) of 
galaxies in the spectroscopic catalog. Keck spectro- 
scopic selection Sps{Isi4:,rd), and combined survey 
selection function Sus{l8iiTfd)- The highest contour 
was normalized to one in each case. The red points 
in the bottom right-hand panel are the galaxies with 
Keck redshifts in the sample. The observed distribu- 
tion of galaxies in the /814 — r^ plane is given by the 
resident galaxy distribution in that plane multiplied 
by Sus{Isi4:: fd)- The solid line in the upper left-hand 
panel is an arbitrary line of constant apparent surface 
brightness. 

Fig. 5. — Disk scale length versus disk rest-frame B- 
band absolute magnitude as a function of redshift for 



galaxies with B/T < 0.5. Filled circles: Galaxies 
with Keck redshifts; Open triangles: Galaxies with 
CFRS redshifts. Long-dashed line: Canonical Free- 
man relation for a disk central surface brightness 
fXOg — 21.65 mag arcsec"^ (Freeman 1970). Error 
bars are 99% confidence intervals. The lowest three 
redshift bins cover nearly the same redshift ranges 
as those in Schade et al. (1996b). The loci of the 
DEEP/GSS galaxies {solid outlines) are used later in 
the coinparison with other samples (Section |3.l| and 
Figure 0) . 

Fig. 6. — Rest-frame selection function Sus{^lB(n 
Rd) {shaded contours) as a function of redshift. The 
highest contour was normalized to one in each red- 
shift bin. Red symbols are galaxies from Figure H 
with Keck redshifts and B/T < 0.5, and the solid 
green line is the canonical Freeman relation. 

Fig. 7. — Loci of different magnitude-size samples: 
DEEP/GSS {solid outline), Schade et al. (1996b) 
{dotted outline), Roche et al. (1998) {dashed outline), 
and Phillips et al. (1997) {dot-dashed outline). The 
DEEP/GSS outlines are the same as those shown in 
Figure |. Note that the DEEP/GSS sample pene- 
trates to lower surface brightnesses at high redshift 
than the other samples. 

Fig. 8. — Raw disk rest-frame B-band central surface 
brightness distributions as a function of redshift un- 
corrected by the survey selection function Sus- Sur- 
face brightnesses were computed using Equation |6[ 
The area under each distribution was normalized to 
one, and the vertical dashed lines mark the value of 
the canonical Freeman disk central surface brightness 
(/iOe = 21.65 mag arcsec"^). The large filled circles 
give the location of the mean surface brightness in 
each redshift bin. The error bar on the location of the 
mean is the 3-sigma standard error. The mean surface 
brightness increases systematically from fiQ^ — 21.7 
mag arcsec"^ at z = 0.2 to hq^ — 20.4 mag arcsec"^ 
at z = 0.80, similar to that seen by Schade et al. and 
Roche et al. 

Fig. 9. — Adjusted disk rest-frame B-band central 
surface brightness distributions as a function of red- 
shift scaled by the survey selection function Sjjs 
following Equation 0. This scaling converts each 
distribution to what it would look like if observed 
with the same selection function as the highest red- 
shift bin. Each galaxy thus contributes the frac- 
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tion Sus{Rdijk I My-fc, 0.90 < z)/Sus{Rd,jk > M,jk,Zi) redshift. 

to these adjusted distributions. Ngai is the number of 
galaxies in each bin for which Sus {Rdi ,, , Mijk ,0.90 < 
z) ^ 0. The vertical dashed lines mark the value of 
the canonical Freeman disk central surface brightness 
(MOb = 21.65 mag arcsec"^). The large filled circles 
give the location of the mean surface brightness in 
each redshift bin. The error bar is the 3-sigma stan- 
dard error on the mean. There is no systematic in- 
crease in the mean surface brightness from z = 0.2 to 
z = 0.8. This shows how important selection effects 
can be in interpreting the magnitude-size relation at 
high redshift. A significant fraction (~ 25%) of galax- 
ies in the highest redshift bin remained close to the 
Freeman canonical surface brightness of /io^ = 21.65 
mag arcsec"^. 

Fig. 10. — Half-light radius and total rest-frame B 
magnitude data for 957 galaxies taken from Burstein 
et al. (1997) and Bender et al. (1992). Up- 
per left-hand panel: ellipticals E {filled circles) and 
dwarf ellipticals dE (open circles). Upper right-hand 
panel: Sa {filled circles), Sab-Sb {open circles) and 
Sbc {pluses). Lower left-hand panel: Sc {filled circles) 
and Scd-Sdm {open circles). Lower right-hand panel: 
Sm {filled circles) and Irr {open circles). The solid 
line is the half light radius-magnitude relation for a 
canonical Freeman disk with a central surface bright- 
ness of 21.65 mag arcsec"^. The Freeman relation is 
provided for the E galaxies only as a reference since, 
strictly speaking, it applies only to disk galaxies. 

Fig. 11. — Half-light radius in kiloparsecs versus total 
rest-frame B-band absolute magnitude as a function 
of redshift for all galaxies with Keck redshifts. Filled 
circles: 0.0 < B/T < 0.5, Open circles: 0.5 < B/T < 
1.0. The solid line is the half-light radius-magnitude 
relation for a canonical Freeman disk with a central 
surface brightness of 21.65 mag arcsec"^. 

Fig. 12. — Direct half-light radius-magnitude compar- 
ison between local BBFN galaxies (dots) and galaxies 
in the DEEP/GSS redshift sample (filled circles). The 
HSB galaxies are shown as open circles. The dashed 
box outlines the magnitude-size locus of the z ~ 3 
galaxies from the sample of Lowenthal et al. (1997), 
and the solid line is the Freeman relation. Compari- 
son with Figure p^ shows that the DEEP/GSS sample 
probably includes different morphological mixtures in 
different redshift bins, ranging from Sm-Irr galaxies at 
low redshift to early-type and HSB galaxies at high 
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Table 1 
Rest-Frame Surface Brightness Distribution Parameters 





No Selection Function 


Adj 


ustment 


With Selection Function 


Adjustment 


z 


^gai 


Mob 




C^(A*Ob) 


Ngai 


A^Ob 


ct(mob) 


0.10 < z < 0.30 


n 


21.7 ± 0.3 




1.2 


5 


20.3 ± 0.6 


1.3 


0.30 < z < 0.50 


18 


21.0 ± 0.3 




1.2 


6 


20.2 ± 0.2 


0.6 


0.50 < z < 0.70 


27 


20.7 ± 0.2 




0.8 


14 


20.3 ± 0.2 


0.8 


0.70 < z < 0.90 


33 


20.4 ± 0.2 




1.1 


31 


20.2 ± 0.2 


1.2 


0.90 < z < 1.10 


41 


19.9 ± 0.2 




1.4 


41 


19.9 ± 0.2 


1.4 
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A. APPENDIX A: Rules for the Spectro- 
scopic Selection Function Sps 

Ideally, Spsihuifd) should just be the ratio of 
the spectroscopic and photometric samples binned in 
both magnitude and size. Unfortunately, this is not 
always the case for two reasons. First, even with large 
bin sizes {Amag =1.0 and A log r^ = 0.5), small- 
number statistics at the bright end cause large fluc- 
tuations in Sps{Iai4:,rd) where it should have been 
flat and close to its maximum value. To alleviate 
this small-number problem, the spectroscopic and 
photometric samples were binned only in magnitude 
(Amag — 1.0). The resulting Sps{Iai4:,rd) selection 
function is thus insensitive to any size-dependent ef- 
fect in the spectroscopic selection. This would be a 
problem if the size distributions of galaxies in the 
spectroscopic sample were markedly different from 
those of the general population as a function of mag- 
nitude. No difference was found in four magnitude 
ranges: 20.5 < hu < 21.5, 21.5 < Isia < 22.5, 
22.5 < /sM < 23.5, and 23.5 < hu < 24.5. Sec- 
ond, the ratio Vts/vl/p of the spectroscopic and pho- 
tometric samples must be calculated in each bin in 
the J814— log rd plane according to the rules shown in 



Table A2. Table A2 describes all the different ways 
in which the selection functions operate on \E'j/ and 
*P to produce *s- "Y" and "N" are used to indi- 
cate whether galaxies in the magnitude-size bin are 
present or not in the ^'s. For simplicity, zeroes and 
ones are used for the selection functions. In reality, 
they can obviously take on any value between and 1 
since they give the probabilities of detecting objects 
photometrically and spectroscopically. 

Cases (3) and (7) are the only two cases where Sps 
is greater than Sup. These cases represent serendipi- 
tous observations of very strong emission-line objects 
which are not detected on broad-band images. Since 
these objects were not included in the spectroscopic 
sample used for this paper. Cases (3) and (7) do not 
apply here. Sps must be carefully interpreted when 
both ^p and '^s sue zero. Take Case (5) for ex- 
ample. Objects in this case are absent from both ^p 
and ^s because they do not exist in Universe and not 
because they could not have been detected in princi- 
ple. It would therefore be wrong to set Sps to zero if 
these selection functions were to be applied to predic- 
tions from theoretical models. The selection functions 
should not erase potential differences between theo- 
retical model predictions and the true galaxy popu- 
lation in the Universe. The best estimate for Sps 



in this case is therefore an upper limit given by the 
value of Sjjp in the same magnitude-size bin. There- 
fore, Sps was set to Sup whenever both '^p and ^s 
were zero. Otherwise, Sps was simply set equal to 

The selection rules discussed so far in this Ap- 
pendix apply to the very general case of a survey 
in which redshifts may not have been obtained for 
all the objects observed and in which redshifts may 
have been obtained serendipitously for objects ab- 
sent from its photometric catalog. The DEEP/GSS 
sample (/8i4 < 23.5) used in this paper was much 
simpler to analyze than such a general survey since 
Keck redshifts were obtained for 100% of the ob- 
jects in the sample, and all objects were present in 
the DEEP/GSS photometric catalog. In this case, 
the spectroscopic selection function Sps in a given 
(/8i4,rd) bin is simply the ratio of the number of 
galaxies observed spectroscopically and the number 
of galaxies in the photometric sample. 
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Table A2 
Spectroscopic Selection Function Sps Rules 



Case 


Universe 


Sup 


Photometric Catalog 


Sps 


Spectroscopic Catalog 




*i7 




*P 




*s 


1 


Y 


1 


Y 


1 


Y 


2 


Y 


1 


Y 





N 


3 


Y 





N 


1 


Y 


4 


Y 





N 





N 


5 


N 


1 


N 


1 


N 


6 


N 


1 


N 





N 


7 


N 





N 


1 


N 


8 


N 





N 





N 
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